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Abstract 


A hybrid model composed of a least square support vector machine (LS-SVM) model and a pressure-incremental model is developed to dispose 
operation conditions of current, temperature, cathode and anode gas pressures, which have major impacts on a proton exchange membrane fuel cell’s 
(PEMFC) performance. The LS-SVM model is built to incorporate current and temperature and a particle swarm optimization (PSO) algorithm 
is used to improve its performance. The optimized LS-SVM model fits the experimental data well, with a mean squared error of 0.0002 and a 
squared correlation coefficient of 99.98%. While a pressure-incremental model with only one empirical coefficient is constructed to for anode and 
cathode pressures with satisfactory results. Combining these two models together makes a powerful hybrid multi-variable model that can predict a 
PEMFC’s voltage under any current, temperature, cathode and anode gas pressure. This black-box hybrid PEMFC model could be a competitive 


solution for system level designs such as simulation, real-time control, online optimization and so on. 
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1. Introduction 


As a clean energy conversion technology, proton exchange 
membrane fuel cells (PEMFC) receive more attention because 
of their low operating temperature, high power density, quick 
start-up capability and long lifetime. PEMFC is an interesting 
technology for the next generation of vehicles, portable units 
and so on [1]. A convenient PEMFC model can help greatly to 
control, simulate, and diagnose its behavior. 

The PEMFC system is a nonlinear, multi-variable electro- 
chemical system that is hard to model. A large number of 
publications on fuel cell modeling [2] target the complicated 
internal phenomena at the molecular level. Among them, two- 
dimensional and more complex three-dimensional, two-phase 
and non-isothermal models have been presented [3-5], these 
had very complicated expressions with some key physical 
parameters that were even immeasurable. These mechanistic 
models usually focused on the electrochemistry, thermodynam- 
ics and fluid mechanics. Typically, they were centered on the 
membrane-electrode assembly (MEA), which could help ana- 
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lyze, design and optimize cell components in the laboratory 
context. However, generally they were not suitable for system 
level research [6-9]. 

An empirical modeling approach is more practical in some 
applications. Researchers can deduce a PEMFC stack’s voltage 
responses without knowing the fuel cell’s complicated inter- 
nal characteristics. An active empirical modeling methodology 
in recent years is based on machine learning theories, such as 
artificial neural networks [10-14] and support vector machines 
(SVM) [15,16]. By mapping the fuel cell voltage as a function 
of various operational conditions, these black-box models agree 
well with experimental data. A common requirement in using 
these modeling approaches is that sufficiently representative data 
should be supplied in the training set to build a multi-variable 
empirical model. However, the numbers of the training data 
needed will increase dramatically when the numbers of input 
variables are increased, and these data may not be available. 
For example, in our previous work [15], 100 experimental data 
points were used to build a two-variable PEMFC voltage model, 
but the number of experimental data points used by Liet al. [16] 
had achieved 1000 to build a five-variable PEMFC tempera- 
ture model. Therefore, although current density, temperature, 
cathode and anode gas pressures are the most important con- 
trollable variables for a PEMFC’s performance, few of those 
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models proposed can deal with all of them simultaneously. One 
important reason is simply the absence of sufficient empirical 
data describing a broad range of operating conditions. 

Herein, a hybrid modeling approach is proposed. In this way, 
only limited numbers of empirical data are needed to build a 
model which can deal with current density, temperature, cath- 
ode and anode gas pressures simultaneously. This model consists 
of two parts: one is an empirical LS-SVM model that concerns 
current density and temperature. Li et al. [16] built a dynamic 
temperature model of a PEMFC by LS-SVM, incorporating tem- 
perature, cooling water speed, hydrogen speed, air speed and 
output power. As a least square version of SVM, LS-SVM inher- 
its the superior generalization performance of SVM. Although 
its precision is slightly lower than SVM, LS-SVM significantly 
reduces computation time. 

The key to obtaining a highly accurate SVM or LS-SVM 
model is to choose a proper set of hyper-parameters, but no 
effective guide lines have been put forward; some recommen- 
dations on how to determine them are quite contradictory and 
confusing. A trial method is used in [15,16], which greatly relies 
on the user’s experience. Therefore, particle swarm optimization 
(PSO) [17] is adopted in this paper to automatically determine 
the best set of hyper-parameters. 

The other part of the hybrid model is a pressure-incremental 
model taking account of both the cathode and anode gas pres- 
sures. This model is developed based on the work of Amphlett 
et al. [18] and Mann et al. [19]. After a slightly simplification, 
the number of empirical coefficients was greatly reduced from 
more than ten to only one. 


2. Theoretical framework 


For a given PEMFC system, the fuel cell terminal voltage 
V is influenced by many operating parameters: cell temperature 
T, cathode oxygen pressure Po,, anode hydrogen pressure Pp,, 
relative humidity a, membrane humidity A, etc. Accordingly, 
fuel cell voltage is given by 


V = fU, T, Poy, Pm, A, a, ...) a) 


Up to now, no model has ever been able to accommodate all 
these operating parameters. With the assumption that channel 
gas is fully saturated and membrane is fully humidified, Eq. (1) 
is simplified as 


V = f(I, T, Po,, Pu) (2) 


Although Eq. (2) is a simplified equation, it is still hard to model 
with traditional means due to its high dimensionality. A way to 
deal with large scale complex systems is to break them into 
independent simple submodels. Thus, we further separate Eq. 
(2) into two parts as 


V = Vp, T) + Vap(APo,, A Pu) (3) 


The first part Vp, denotes a LS-SVM model, which predicts cell 
voltages at different currents and temperatures under a constant 
cathode gas pressure P8, and a constant anode gas pressure PÈ 5 
We refer to these constant pressures as reference pressures. The 
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Fig. 1. Framework of the hybrid pressure-incremental LS-SVM model. 


second part Vap denotes a pressure-incremental model, which 
predicts voltage increment caused by oxygen pressure increment 
A Po, and hydrogen pressure increment A Py,. The structure of 
our proposed hybrid pressure-incremental LS-SVM model is 
illustrated in Fig. 1. 


3. Optimized LS-SVM model 


A support vector machine is a novel and powerful tool based 
on statistical learning theories. It was originally developed at 
AT&T Bell Laboratories by Vapnik [20] for classification in 
various domains of pattern recognition, then expanded success- 
fully to deal with regression problems more recently. The SVM 
model possesses a high degree of precision and does not require 
a pre-knowledge of the fuel cell. Comparing to artificial neural 
networks, the SVM model has a superior capability of general- 
ization and it is also more robust. LS-S VM proposed by Suykens 
and Vandewalle [21] is a least square version of the standard 
SVM. Compared to SVM, LS-SVM significantly reduces the 
computation time with a tiny precision loss. Taking advantage 
of the high computing efficiency of LS-SVM, PSO strategies are 
used to automatically optimize performance. 


3.1. Modeling a PEMFC by LS-SVM 


Building a LS-SVM model need three steps: 


e Preparing training data. 

e Selecting optimal LS-SVM parameters to train the LS-SVM 
model. 

e Predicting with the LS-SVM model. 


In this study, experimental data provided by Laurencelle et al. 
[22] was used to generate the training data for the LS-SVM 
model. This fuel cell was composed of 36 cells; each cell with a 
232 cm? active area, graphite electrodes, and a Dow membrane. 
Air pressure and hydrogen pressure were both regulated to 3 atm. 
Training data were obtained at 24°C, 31°C, 39°C, 56°C and 
72°C and current densities from 0 mA cm? to 1000 mA cm7?. 
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A difficult and vital step in obtaining a highly accurate LS- 
SVM model is to choose a proper set of hyper-parameters y and 
o. y determines the trade-off between model errors and model 
complexity, and o is a internal kernel parameter. However, a few 
recommendations in the literature on how to determine them 
are quite contradictory and confusing. Past literature [15,16] 
used a try and trail approach to determine the optimal hyper- 
parameters, which needed much experience and patience by 
developers. Therefore, a PSO algorithm was adopted to auto- 
matically search for the best set of hyper-parameters. 


3.2. Optimizing the LS-SVM model by PSO 


PSO is an evolutionary computation technique proposed by 
Kennedy and Eberhart [17]. PSO is an extremely simple algo- 
rithm that requires only primitive and simple mathematical 
operators, and is computationally inexpensive in terms of both 
memory requirements and time. PSO is an algorithm based 
on population (“swarm”). Each individual, called a particle, 
flies in the D-dimensional problem space with a velocity that 
is dynamically adjusted according to the flying experience of 
the individual and its colleagues. Standard PSO algorithm is 
described as follows: 


v; = wv; + c;Rand()(p; — xi) + e2rand()(p, — xi) (4) 
X; = X; + Vi (5) 


where x; represents the current position of particle i, v; is its 
velocity, p; is the best position it experienced. p, represents 
the best position among all particles in the swarm. w is iner- 
tia weight, which regulates the trade-off between the global 
exploration and local exploitation abilities of the swarm. The 
acceleration constants c; and c2 represent the weight of the 
stochastic acceleration terms that pull each particle toward self 
best position p; and swarm best position pg. Rand() and rand() 
are two random functions with range [0,1]. The basic idea of 
PSO is that every particle’s action v; is excited by self cognition 
p; and social experience pg to fly to the globe best position. 

Using model’s mean square error (MSE) as fitness, the opti- 
mization problem is defined as follows 


N 
; 1 a2 
minJ(y, o) = MSE = X > Vi — Vi) (6) 


where N is the numbers of experimental data, V; is experimental 
data, and V; is predicted data. The process of implementing PSO 
algorithm to solve this optimization problem is as follows: 


1. Initialize: initialize a swarm of particles with random posi- 
tions and velocities in the three-dimensional (y, o). Number 
of particles is 10; w is set to 0.4; cı and cp are all set to 2. 
Maximal generations number is 50. 

2. Evaluate each particle’s fitness: train the LS-SVM model by 
current position (y, o);, and then evaluate its fitness by Eq. 
(6). 

3. Update each particles’ best fitness value p; according to the 
result of step 2. If p; is best than the swarm’s best fitness 
value pg, then reset pg to pj. 


0 —_. ; 
500 ~~~ n _ 40 
>< 


Current density (mA/cm2) 1000 20 temperature (°C) 


Fig. 2. Voltage predicted by LS-SVM model. Asterisk data is the experimental 
data. 


4. Change the velocity and position of the particle according to 
Eqs. (4) and (5), respectively. 
5. Loop to step 2 until maximal generations number is reached. 


Above program is written in matlab language and implemented 
on Matlab 7.0. The minimum MSE of 0.0002 is found at 


(y, ©) = (650, 1.38) 


As shown in Fig. 2, the LS-SVM model with optimal hyper- 
parameters fits the experimental data quit well. The mean 
squared error is 0.0002 and the squared correlation coefficient 
is 99.98%. 

To make the model incorporate cathode and anode pressures 
besides current and temperature, a big problem should be con- 
sidered firstly: if increasing the input variables form (/, 7) to 
(I, T, Po,, Py,) to build a high dimensionality LS-SVM model, 
a huge amount of training data would be needed. Assuming n 
training data is needed for one variable, then n? training data is 
needed to build a 2-variable model, and n* training data is needed 
to build a 4-variable model. The numbers of training data per 
variable increase exponentially with the number of variables to 
maintain a given level of accuracy, and these data may not be 
available. To solve this problem, a pressure-incremental model 
was developed to expand our proposed LS-SVM model to other 
pressures. 


4. Hybrid pressure-incremental LS-SVM model 
4.1. Pressure-incremental model background 


An influential one-dimensional PEMFC model was proposed 
by Amphlett et al. [18], and Mann et al. [19] further developed it 
to a generalized form. Herein, a simplified pressure-incremental 
model is developed based on their work. When current is drawn 
from a fuel cell, the cell voltage Vee is decreased from its 
equilibrium thermodynamic potential Enernst (open circuit volt- 
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age). This voltage drop includes activation loss act, ohmic loss 
Nohmic and concentration loss neon. The concentration loss is 
only notable when the fuel cell works with an extremely high 
current [23]; hence it can be neglected under normal conditions. 
The basic expression for the cell voltage is: 


Veet! = Enernst — Nact — Nohmic (7) 


In Section 4.2, we firstly demonstrate the relationship between 
pressure and each of the terms in Eq. (7), and then describe how 
the pressure change affects them. 

The first term is reversible thermodynamic potential Enerst, 
which is described by the Nernst equation. Using literature val- 
ues for the standard-state entropy change, the expression is: 


Enemst = 1.229 — (8.5 x 1074) - (T — 298.15) 
+(4.308 x 1075) - T - (In Py, + 0.5 -In Po,) (8) 


where T is cell temperature (K), Py, the partial pressure of 
hydrogen at the anode catalyst/gas interface (atm), and Po, the 
partial pressure of oxygen at the cathode catalyst/gas interface 
(atm). 

The second term is activation overvoltage Nact, which is com- 
posed of anode and cathode activation overvoltages. Grouping 
them together, we can get a single expression as: 


Nact = E1 + oT + &T In Co, + &4T ln T (9) 


where Co, is the concentration of dissolved oxygen at the 
gas/liquid interface, J is the cell current (A), and those four &; 
(i= 1-4) are coefficients. 

The third term in Eq. (7) is ohmic polarization, which is 
caused by resistances in electron and proton transfers and other 
contact resistances. It can be expressed as: 


Nohmic = IR™™! = I(E5 + £T + 71) (10) 


where &; (i= 5-7) are coefficients. It can be concluded from Eq. 
(10) that the ohmic overvoltage is independent of the pressure, 
which has no influence on our pressure-incremental model and 
could be neglected. Afterwards, a pressure-incremental model 
with only one coefficient is developed. 


4.2. Pressure-incremental model 


As shown in above discussion, hydrogen and oxygen pres- 
sures influence the cell voltage through the terms Py,, Po, and 
Co,. The partial pressure of hydrogen at the gas/liquid catalyst 
layer interface is [18] 


Pp, = Panode + [1 — 0.5xo] (11) 


The partial pressure of oxygen at the gas/liquid catalyst layer 
interface is [18] 


0.291 -1I 
Po, = Peathode «| 1 — afio — 2h - exp (am )| 
(12) 


Fig. 3 demonstrates the values of the exponent term 
exp(0.291-1/T°837) at 40°C and 80°C. From the figure it can 
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Fig. 3. Values of the exponent term at 40°C and 80°C. 


be indicated that the difference with and without this term is no 
more than 0.25%, which is tiny enough to be neglected. 

Hence, the partial pressure at the gas/liquid catalyst layer 
interface is simplified as 


Po, = Peathode * [1 — xfitg — xR] (13) 
The third term is the oxygen concentration Co, [18]: 

P 
Co; = (14) 


2 ™ (5.08 x 10% - exp(—498/T) 


Let PS, and PO, standing for the given reference partial pres- 
sures of oxygen and hydrogen, respectively, under which the 
LS-SVM model predicts the fuel cell voltages at different cur- 
rents and temperatures. The oxygen increment and the hydrogen 


pressure increment are defined as the proportional form: 
0 0 
Po, = ko, Po,» Pu = km Ph, (15) 


where ko, and kp, are ratios that represent the pressure incre- 
ments to PS, and PË,- With the substitution of Eq. (15) into Eq. 
(8), the thermodynamic potential can be rewritten as: 


Enernst = 1.229 — (8.5 x 1074) - (T — 298.15) 
+(4.308 x 107°). T - (In ke, Py, + 0.5 - In ko, P8,) 


Fherg + (4.308 x 1075) -T - (n kp, + 0.5 - In ko,) 


with F etisi 
= 1.229 — (8.5 x 1074) - (T — 298.15) 
+(4.308 x 107°). T - (In PẸ, + 0.5- In P8) (16) 


where E? mnst is the thermodynamic potential under the reference 


pressures Po, and Puy: With the substitution of Eqs. (15) and 
(14) into Eq. (9), the total overvoltage can be rewritten as 


Nact = &1 + &T + &T Inko, Co, + ¿4T In] 
= £1 + T+ £3T In C8, + &4T In I + &T Inko, 
=n? +é&Tlnko, withn?, 
= $1 + éT + &3T In CQ, + &T inl (17) 


where n? is the overvoltage under the known pressure P°. With 
the substitution of Eqs. (16) and (17) into Eq. (7), the fuel cell 
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Table 1 

Errors between actual voltage and predicted voltage 

No. Po, /ko, (atm) Pp, / ky, (atm) Actual voltage (V) Predicted voltage (V) Error (%) 
1 1.40/— 2.50/— 0.824 - - 

2 3.10/2.2143 3.10/1.2400 0.851 0.8586 0.92 

3 0.60/0.4286 2.00/0.8000 0.781 0.7924 1.38 

4 0.60/0.4286 3.10/1.2400 0.792 0.7988 0.83 

5 3.10/2.2143 2.00/0.8000 0.847 0.8473 0.04 


voltage can be rewritten in an incremental format as: 
V = Vp + Vap with Vap 
= §T In ko, + (4.308 x 1075). T - (in ky, + 0.5 - In ko,) 
(18) 


where Vp, is the voltage predicted by LS-SVM model under 
the reference pressures Pe, and Puy: Eq. (18) is our proposed 
pressure-incremental model. By neglecting exponent term in Eq. 
(12), above algebraic manipulations greatly simplifies the elec- 
trochemical model by reducing the empirical coefficients £1, £2, 
&3, &4, &5, E6, &7 and several other numeric empirical constants 
to only one coefficient 3. All those non-pressure-dependent 
terms are effectively aggregated into a single measured refer- 
ence pressure. &3 is a coefficient which relates only to chemical 
parameters of the cathode reaction and does not varies with 
operation conditions. Simple experiments can easily determine 
it. Mann et al. [19] provided a range of (7.6+0.2) x 1075, 
and the middle value of 7.6 x 1075 is adopted in this 
paper. 

To verify our proposed model, experimental data proposed 
by Amphlett et al. [24] is used, as listed in Table 1. These data 
are obtained with temperature of 70°C and current density of 
1.03 Acm~. The first group is used as reference data for our 
predictions, and the other four groups as data to verify our pre- 
dictions. Based on the reference data, Eq. (14) can help predict 
voltages under other pressures. Although our model is a sim- 
plified form of Amphlett’s model, the modeling results are still 
highly consistent with the same measurements. 

Fig. 4 shows the fuel cell voltages at 70°C and 6.66 A. It is 
shown that the oxygen partial pressure bears a stronger relation- 
ship to the fuel cell voltage than the hydrogen partial pressure 
does. Thus, it will improve a fuel cell’s output power more effi- 
ciently if the cathode pressure is increased rather than the anode 
pressure. 


4.3. Hybrid model 


In Section 3.2, the optimized LS-SVM model has been 
built to predict the voltages at air pressure and hydrogen pres- 
sure are both 3 atm with temperature from 24°C to 72°C. 
Using this prediction as reference data, our pressure-incremental 
model can be made to take the prediction at other cathode 
and anode pressures. As an example, the predictions at cath- 
ode/anode pressure at 1/1 atm and 5/5 atm are demonstrated in 
Fig. 5. 


Channel partial pressure (atm) 


Fig. 4. Fuel cell voltages at 70°C and 6.66 A. The dot line represents the refer- 
ence cell voltage under the working conditions in Table 1. In the upper curve, 
Po, varies from 1 atm to 6 atm with Py, kept at 2.5 atm. In the lower curve, Py, 
varies from 1 atm to 6 atm with Po, kept at 1.4 atm. 


In order to get exact prediction results for different PEM- 
FCs, £3 should be calibrated based on corresponding benchmark 
experimental data. Since &3 is a constant for a given fuel cell, only 
small numbers of experimental data should do this procedure 
well. 
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Fig. 5. Voltage predicted by the hybrid model under different pressures with 
temperature from 25°C to 70°C. 
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5. Conclusions 


In this study, we have developed a system-level hybrid model 
of a PEMFC incorporating current, temperature, cathode and 
anode pressure. This model consists of two parts. One is an 
empirical LS-SVM model that concerns current and temper- 
ature. After optimization by a PSO algorithm, this LS-SVM 
model possesses excellent agreement between the predicted 
and experimental results. Off-line optimizing of the LS-SVM 
model took 68 min on a Pentium IV 1.6GHz computer with 
512 MB RAM, but the predicting time was no more than 1 ms. 
The other one is a pressure-incremental model that concerns 
cathode and anode pressures with only one empirical parame- 
ter. With little increase of benchmark experimental data in the 
two-variable LS-SVM model, the hybrid model could deal with 
all four operational variables. This hybrid black-box model- 
ing methodology is most applicable for system level research 
such as simulation, real-time control, online optimization and 
so on. 

One flaw of our hybrid model is that it should be used under 
modest changes in pressure, as significant changes in pressure 
will affect both the voltage and current, thus moving the fuel cell 
away from the reference conditions under which the simplified 
pressure-model was calibrated. It also will be invalid in low 
humidity or under extremely high current density operations. 
Hence, further studies are still needed. 
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